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Abstract 

We reconsider Archimedes evaluations of several square roots in 
'Measurement of a Circle'. We show that several methods proposed 
over the last century or so for his evaluations fail one or more criteria 
of plausibility. We also provide internal evidence that he probably used 
an interpolation technique. The conclusions are relevant to the precise 
calculations by which he obtained upper and lower bounds on ir. 
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1 Introduction 

We do not know how Archimedes proved his famous result 

265 r 1351 
153 <V3< W 

mentioned in 'Measurement of a Circle', but a number of possible methods 
have been suggested. They may be found on the web and elsewhere; see [TJ HI 
El EH] • A good candidate should satisfy the following criteria. 

• It should only use methods that were known at the time of Archimedes; 

• It should be possible to replace 3 by any other positive integer; 

• It should be as short and elementary as possible; 

• It should not involve extraordinary ingenuity or tricks. 
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Some of the proposals discussed by Heath in [H pp. 51, 52] fail the first test, 
unless some much later authors were describing methods that had been known 
centuries before their own time. 

Section 4 describes a new 'interpolation' method that meets all four criteria, 
and that, if correct, would further confirm Archimedes' understanding of limits. 
In Section 6 we prove that Archimedes' lower bound cannot be obtained by 
Hero's method. Section [7] discusses a number of much harder square roots that 
Archimedes calculates in the course of obtaining his even more famous upper 
and lower bounds 

3±f < 7T < 3f. 

Several quite difficult new issues emerge, but it is again likely that Archimedes 
used an interpolation method to calculate the square roots. 

Before starting we mention three standard problems in this type of enterprise. 
Euclid's Elements survives in a variety of substantially different versions, all 
of which are copies made hundreds of years after he died. Although he must 
have written it around 300 BC, little is known about him personally. For this 
reason the word Euclid should be regarded as a label rather than a name. See 
[3j Chap. 2, 4] and [6, pp. 354-361] for more detailed accounts of this issue. 
Two comprehensive sources of the Elements are [8] and [5], but a valuable 
summary may also be found at 0. 

Archimedes of Syracuse (c. 287-212 BC) is a less shadowy figure than Eu- 
clid, but our knowledge of his mathematical work is similarly indirect. His 
'Measurement of a Circle' survives via numerous copies and translations of a 
ninth century document written in Constantinople, [3 pp. 25-27]. It is also a 
part of the famous Archimedes palimpsest, transcribed in the tenth century, 
overwritten in the thirteenth century and discovered in 1906 by Heiberg. The 
palimpsest is now being studied by means of sophisticated imaging techniques 
at the Walters Art Museum in Baltimore, [2]. 

The second problem is that with the knowledge of later developments, one 
can easily read more into an ancient text than is there. In particular, it is 
tempting to assume that an author was aware of what now seems to be an 
obvious development of an idea in his surviving work. It has, for example, 
been argued that Euclid and his immediate successors did not understand the 
fundamental theorem of arithmetic in the way that we now do, [12] . One must 
constantly be aware of the dangers of falling into such a trap. 

The third problem is finding a satisfactory compromise between notation and 
concepts that Archimedes would have understood and our very different way of 
expressing the problems of interest. This article uses modern terminology for 
the calculations, including the free use of negative numbers, but only in circum- 
stances in which there seems to be no problem with rewriting the arguments 
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appropriately. 

Throughout the paper Roman letters are used to denote natural numbers and 
Greek letters to denote fractions or other real numbers. All numerical values 
are rounded down so that the digits presented are digits of the exact value. 

2 Brute force 

The brute force method involves calculating m 2 and 3n 2 for all m, n up to 
1351, and then finding m, n such that m 2 — 3n 2 is very small. If one calculates 
one square per minute, producing such a table would take twenty-four hours, 
and the table would then be available for other uses. The squares need not be 
calculated ab initio because of the identity (n + l) 2 = n 2 + 2n + 1; this has 
the independent advantage that one can detect occasional errors by examining 
the result for every n that is a multiple of 10. Napier and others completed 
much more arduous tasks before computers were invented. The fact that the 
method fails the third of our tests is not decisive, in spite of the complexities 
of Greek arithmetic when dealing with numbers much bigger than a thousand. 
Archimedes might have used this method, but it would have been much more 
laborious than the interpolation method that we describe in Section HI 

Table Q] shows a fragment of a systematic calculation of the squares of the first 
thousand natural numbers. 
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X 
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= 44652 
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123 2 


= 45387 




+247 








+741 


124 2 = 


15376 


3 


X 


124 2 


= 46128 



Table 1: Calculations of squares 

Table [2] presents a selection of the results of the calculations indicated above. 
In the second column the upper bounds to a/3 are displaced slightly to the 
right. 



3 



5 2 - 3 x 3 2 = 
7 2 - 3 x 4 2 
19 2 -3 x ll 2 = 

26 2 - 3 x 15 2 
71 2 -3 x 41 2 = 

97 2 - 3 x 56 2 



2, 

1, 
2, 
1, 
2, 
1. 
2, 
1. 
2. 
1. 



5/3 = 1.6, 
7/4 = 1.75, 



71/41 ~ 1.731707, 
97/56 ~ 1.732142, 



19/11 ~ 1.727272, 
26/15 ~ 1.733333, 



265 2 - 3 x 153 2 = 



265/153 ~ 1.732026, 
362/209 ~ 1.732057, 



362 2 - 3 x 209 2 



989 2 - 3 x 571 2 = 



989/571 ~ 1.732049, 
1351/780 ~ 1.732051. 



1351 2 - 3 x 780 2 



Table 2: Results of brute force calculations 



3 Pell's equation 

It is often suggested that Archimedes might have solved the equations a 2 — 3b 2 = 
m for m = 1 and m = — 2 systematically; it should be mentioned here that 
a 2 — 3b 2 = — 1 has no integer solutions. We argue that this is not plausible on 
historical grounds. 

The equation a 2 — cb 2 = m is often called Pell's equation, particularly when 
m = l, because of an unfortunate error; it has nothing to do with Pell. It is 
obviously equivalent to (a/b) 2 — c = m/b 2 , so any integer solution with b large 
and m small yields an accurate rational approximation to y/c. This method 
is indeed one of the best ways of approximating y/c, and the main problem is 
finding evidence that Archimedes would have been able to solve the equation 
when c = 3. The situation is as follows. 

It is known from cuneiform clay tablets that the Babylonians were able to 
calculate square roots accurately long before the time of Archimedes, but there 
is no contemporary evidence of how they did this. In the fifth century AD 
Proclus claimed that the Pythagoreans calculated the square root of 2 by using 
the following fact, which I have rewritten in modern notation. If a n , b n satisfy 
a 2 — 2bf l = ±1 and one puts a n+ i = a n + 2b n , b n+ i = a n + b n then a 2 +1 — 26 2 +1 = 
=Fl. Starting from a\ — b\ — 1 this generates a sequence of steadily improving 
upper and lower bounds alternately. This is, in fact, a systematic procedure 
for solving Pell's equation when c = 2 and m = ±1. 

Proclus's claim may have been based on guesswork or on documentary sources 
that are now lost. One has to be careful not to accept it uncritically, because 
he would have known the much later work on such equations by Diophantus of 
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Alexandria, written in the third century AD. However he quotes the identity 

(a + 2b) 2 -2(a + b) 2 = -(a 2 - 2b 2 ) 

from the Elements, Book 2, Prop. 10. This makes a reasonable case for be- 
lieving Proclus, but it does not imply that Archimedes was able to generalize 
the procedure to c = 3. Indeed there is no evidence to support this. Solving 
such equations might now seem elementary, but the evidence is that it took 
many hundreds of years before the relevant theory became routine knowledge; 
see P21 Chap. 5B]. Some of the identities needed look very artificial if one does 
not express them in term of quadratic arithmetic in the ring Z(a/c), within 
which the equation can be rewritten in the form (a + by/c)(a — byfc) = m. 



4 An interpolation method 

We describe the method in modern notation for the convenience of the reader, 
and explain its relationship with the use of continued fractions in the next 
section. The method differs from others in that it starts with two approxima- 
tions rather than one, and passes by an elementary calculation to a new pair 
of approximations. 

The only result that we use is the fact that 

v y . v v + y y 

— < — implies — < < — . (1) 

U X ' U U + X X 

We have not found this in the Elements, but it can be deduced from Book 5, 
Prop. 12, whose proof is logical rather than geometric. It may also be proved 
by modifying the proof of Prop. 12 appropriately. 

One can obtain an insight into ([I]) by considering the slope of the diagonal 
of a parallelogram that has two adjacent edges with slopes v/u and y/x as 
in Figure [TJ It is not obvious whether one can convert this use of coordinate 
geometry into a form that Euclid or Archimedes could have understood. 

An alternative proof of (PQ) involves rewriting the hypothesis in the form v x < 
yu and interpreting this as comparing the areas of two rectangles. The first 
inequality in the conclusion is rewritten in the form v(u + x) < (v + y)u and 
is also interpreted in terms of areas. The result is then obvious by inspecting 
Figure |2j The second inequality has a similar proof. This type of transforma- 
tion is to be found in the Elements, Book 6, Prop. 14, and ([I]) can be proved 
from Prop. 14. 

Given two fractions a = v/u and (3 = y/x such that a < < (3, one puts 

v + y 

7 = —— ■ 

u + x 
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Figure 1: Interpolation using slopes 



yu 




vu 


vx 



Figure 2: Interpolation using areas 

One then replaces a by 7 if 7 < v^3, or equivalently if (v + y) 2 < 3(u + x) 2 , 
and replaces (3 by 7 if 7 > \/3, or equivalently if (v + y) 2 > 3(u + x) 2 , to obtain 
a new and better enclosure of \/3. This is repeated until the desired accuracy 
is achieved. 

The basic idea of this method, obtaining steadily better upper and lower 
bounds on a number by a systematic algorithm, was used by Archimedes on 
many occasions. See, for example, Section [71 

A list of the first sixteen iterations of the interpolation method starting from 
the pair 1, 2 is shown in Table [3] It may be seen that Archimedes' upper and 
lower bounds are obtained by this method. However, his lower bound is by no 
means the best in the list. There is not enough information to explain this, 
but the following is plausible. He might initially have stopped calculating at 

265 r 362 

< V3 < . 

153 209 
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< 




7 
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19 
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45 
26 


< 




26 
15 


71 
11 


< 




26 
15 


71 
il 


< 


V3< 


97 
56 


168 
97 
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y/3< 


97 
56 


265 
153 


< 




97 
56 


265 
153 


< 


y/3< 


362 
209 


627 
362 


< 


V3< 


362 
209 


989 
571 


< 


Vs< 


362 
209 



989 ^ AT ^ 1351 
571 V 780 ' 



Table 3: The interpolation method 

If so, he would later have discovered that the upper bound on y/3 was not strong 
enough to obtain the lower bound that he desired on tt and then proceeded 
further with the iteration above, but only recorded the improved upper bound 
to a/3, because he did not want to repeat more calculations than necessary. 

The interpolation method is systematic and completely elementary, apart from 
depending on being able to determine whether a fraction a/6 is greater than 
or less than y/3. This is equivalent to whether a 2 — 3b 2 is positive or negative. 
This can be decided by evaluating the squares involved, which is quite labo- 
rious when a and b are large. However, mathematicians of all eras up to the 
invention of pocket calculators were very skilled when carrying out arithmetic 
calculations. 

Dominus has described a method that yields the same upper and lower bounds 
as those obtained by the interpolation method, [I]. It depends on producing 
a list of good approximations to y/3 with small numerators and denominators 
and then finding a pattern that can be extrapolated. It is finally necessary to 
check that the extrapolated values are indeed good approximations. The main 
problem with this method is that if one applies it to more complex problems, 
such as the determination of the square roots of larger numbers, no obvious 
pattern is likely to emerge. On the other hand the interpolation method always 
works because of its obvious similarity to the interval bisection proof of the 
intermediate value theorem. 



5 Continued fractions 

The fraction a/b can only be close to y/3 if m = a 2 — 3b 2 is small relative to 
b 2 . This proves that it should be possible to find a connection between any 
method of approximating y/3 by fractions and solutions of Pell's equation or 
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some quadratic Diophantine generalization. Pell's equation is known to be 
connected in turn to the theory of continued fractions, [2]. In this section we 
show that the method of continued fractions is a special case of the interpolation 
method. This by no means implies that the theory of continued fractions 
underlies the calculations of the last section, which are elementary and rigorous 
in their own right. There is no evidence that the use of continued fraction 
expansions was known in classical Greece, in spite of the fact that determining 
the highest common factor of two numbers by using the Euclidean algorithm 
can be interpreted that way. 

Instead of providing a complete analysis, we consider the particular case in 
which some positive real number x is approximated by 

1 

x ~ ai H ; — 

where a r are all positive integers. Some elementary algebra yields 

1 a + /3a 4 

a i ^ 1 — - ; — (2) 

a 2 H fx 7 + oa± 

where a, (3, 7, 5 are certain polynomials in ax, a 2 , 03. Moreover 

— — ax H (3) 

7 a 2 

and 

f = ax + — ^r- (4) 

The validity of ([2]) follows by considering the map 

1 

04 — > ai H i — 

a 2 + 



an- 



as the composition of a sequence of fractional linear transformations on R U 
{00}. The equations ([3]) and (jlj) can then be obtained by putting 04 = 0, 00 
respectively. These formulae show that if the positive integers a 1; a 2 , 0.3 have 
already been chosen by the continued fraction method to approximate x as 
well as possible, then the choice of the next positive integer 04 is of exactly the 
type discussed in Section HI 



6 Hero's method 

In the first century AD Hero of Alexandria described a method which converges 
more rapidly than the interpolation method. Perhaps over-generously, we will 
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use the term 'Hero's method' to refer to any algorithm for calculating square 
roots that relies upon the fact that the arithmetic mean of two numbers is 
greater than their geometric mean, which in turn is larger than their harmonic 
mean. The AM-GM inequality is mathematically equivalent to the GM-HM 
inequality and is a special case of the Elements, Book 5, Prop. 25. It may be 
proved by observing that one can put four disjoint rectangles of size a x b inside 
a square with edge length a + b, as in Figure [3J These means were known to 
the Pythagoreans and were mentioned in Plato's Timaeus 35, 36. 



Figure 3: The AM-GM inequality 

Hero's original method may be found in his Metrica, which was lost for many 
centuries, but eventually rediscovered in 1896 in Constantinople by R Schone. 
An account of Hero's mathematical work may be found in [3 Chap. 18] and 
other important contributions of his are described in [10]. Applying Hero's 
method to the numbers a and b = 3/a, whose geometric mean is a/3, one 
obtains 

, a + 3/a r- 
a! = — ^— > V3. 

Starting from a = 5/3 this can be applied twice to yield 

26 i— 1351 r 

— > V3 and then — — > V3. 
15 780 

In spite of this, it is not possible to use arithmetic and harmonic means to 
obtain Archimedes' lower bound |||, because of the following theorem; note 
that 3 is not a factor of 265. This proves that Archimedes did not use Hero's 
method to obtain his lower bound. 
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Theorem 1 Let x, y be two distinct positive rational numbers with xy = 3 
and let a, h be their arithmetic and harmonic means respectively. Then ah = 3 
and h < \/3 < a. Moreover 3 must be a factor of the numerator of h. 



Proof The formulae 



x + 3/x 
a = , h 



l/x + x/3 



yield ah = 3. The second statement of the theorem follows from the fact that 
the arithmetic mean of two numbers is greater than their geometric mean, 
which is 3 in this case, while their harmonic mean is smaller. 

In order to prove the final statement we put x = n/m where m, n are relatively 
prime, and consider the resulting expression 

Gran 



3m 2 + n 2 ' 



in which the numerator and denominator may have some common factors. One 
needs to consider the effect of removing these. If 3 is not a factor of n then 
it is not a factor of the denominator, so the 3 is the numerator cannot cancel 
out. If 3 is a factor of n then it cannot be a factor of m, so 3 is a factor of the 
denominator but 3 2 is not. On the other hand 3 2 is a factor of the numerator, 
so after reducing h to its lowest terms, the numerator must retain a factor 3. 
□ 

In his account of Hero's work, Heath describes a complicated mixed strategy 
for obtaining Archimedes' lower bound, [TJ pp. 324, 325]. He first applies 
Hero's method starting from a = 1 to obtain 

97 R 168 m 

56 <VI< lrT (5) 

in three steps. Further iterations of Hero's method from this starting point do 
not yield Archimedes' upper or lower bounds. Then he observes that if one 
applies the interpolation method to the two fractions in (jSJ) one obtains 

97+ 168 _ 265 
56 + 97 ~ 153 

and leaves the reader to infer that this is a lower bound. Finally he re-applies 
Hero's method starting from a = 5/3 to obtain Archimedes' upper bound. 
Although he does eventually obtain the result required, he does so by using 
Hero's method with two different starting points and provides no reason why 
one should then apply the interpolation method to the particular pair of frac- 
tions in §5$) rather than to some other pair. It therefore fails to satisfy our 
third and fourth criteria for a satisfactory solution to the problem. 
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7 Upper and lower bounds on n 



Archimedes' upper and lower bounds 



3^ < 7T < 3i 



on 7i have been discussed in some detail by Heath, [7, pp. 50-56]. We follow his 
account of Archimedes' text, but draw significantly different conclusions about 
the often unstated numerical methods that Archimedes used at various stages 
in the calculation. 

Archimedes' upper and lower bounds were obtained by approximating the edge 
lengths of regular polygons that have 6, 12, 24, 48, and finally 96 edges. There 
are three reasons for his residual error: he started with upper and lower bounds 
to \/3 rather than the exact value, at each stage he obtained upper or lower 
bounds on certain further square roots, and he stopped at the 96-gon stage. 
Numerically his final result is 



With current computational facilities, one can assess the effect of the inexact 
values of Archimedes' square roots as follows. The exact bounds obtained by 
calculating the total perimeter of the 96-gon are 



which translate into the numerical bounds 

3.141031 < 7T < 3.142714. 

By Taylor's theorem, each further doubling of the number of sides reduces both 
errors by a factor of 4. At each stage the error in the upper bound is about 
twice as big as the error in the lower bound - until the rounding errors become 
important. 

Archimedes relied on the following theorem to pass from one polygon to another 
with twice the number of sides. 

Theorem 2 (Elements, Book 1, Props. 5 and 32) Let 9 n be the internal 
angle of the right-angled triangle with edge lengths a n , c n and b n = \/ a\ + c\, 
as shown in Figure^ Then the triangle with edge lengths 



3.140845 < 7T < 3.142857. 



96sin(7r/96) < vr < 96 tan(7r/96), 




has internal angle 9, 



n+l — 



n /2. 
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Figure 4: Halving the internal angle 



Applying Theorem [2] repeatedly, Archimedes had to evaluate seven square 
roots, in addition to y/3, but did not explain how he did this, [TJ pp. 55, 56]. 
Five of the seven square roots are approximated, above or below as necessary, 
by an integer plus a fraction whose denominator is 4 or 8. Any explanation 
of how he performed the calculations must take this fact into account. It is 
probably relevant that dyadic fractions had already been used in measurement 
by the ancient Egyptians. One of his bounds was 

7 = V571 2 + 153 2 > 591|. 

A numerical computation of the LHS yields 7 ~ 591.14296, which confirms his 
result and indicates its level of accuracy. 

It is likely that Archimedes proceeded as follows. First he found the integer part 
n of 7. This could have been done by using Pythagoras' theorem graphically, 
i.e. by drawing the right angled triangle whose two shorter edges had lengths 
571 and 153. This produces the approximate value 590 for the hypotenuse 
immediately, after which a few calculations yield 591 < 7 < 592. He then 
applied the following interpolation method to obtain more accurate upper and 
lower bounds. He bisected the interval [591,592] and chose the subinterval 
containing 7 by calculating the square of the midpoint; after repeating this 
procedure three times, he chose the left hand end of the resulting subinterval 
because he needed a lower bound to 7. If this is how he proceeded, he would 
have obtained the results in Table HJ which yields (59l|) 2 < 571 2 + 153 2 , or 
equivalently 591 1 < 7. 

The squares in the above interpolation could have been calculated directly, but 
starting from the fourth line Archimedes might have used the formula 

which would have reduced his labour considerably, because at each stage a 2 , b 2 
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7 2 = 571 2 + 153 2 = 349450 



590 2 = 348100, 
592 2 = 350464, 
591 2 = 349281, 



(591|) 2 = 349872± 



■i • 



(591±) 2 = 349576^, 
(59l|) 2 = 349428§§. 



Table 4: Calculations to determine 7 

and a — b are all known. The formula (E]) is contained in the Elements, Book 
2, Prop. 9 or 10. 

We next discuss the first of the two exceptional cases, in which Archimedes 
recorded the inequality 



a 



V1823 2 + 240 2 < 1838^. 



A numerical computation yields a ~ 1838.7302, so according to the above 
method Archimedes should have recorded the better upper bound 1838|. It 
is plausible that he obtained this value by the interpolation method described 
above, and then adjusted it at the next step, where he used the simplification 

1823 + 1838^ _ 3661^ _ 1007 
240 ~ 240 ~ 66 

This sufficed for his purposes and would have reduced the labour involved in 
later calculations. Unfortunately the fraction 

3661 1 
240 

cannot be simplified. 

The second exceptional case must have a different explanation. Archimedes 
obtained 



P = V1007 2 + 66 2 < 1009| (7) 
which led him after one further calculation to 

96 x 66 ^ (^ 3.140845). (8) 

20171 71 v ) v ) 
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It is suggested that he originally proved that (3 < 1009| by the interpolation 
method, and then continued to obtain 

96 x 66 Q 
tt> ^>3|j = 3.140625. 

Not satisfied with this, he modified the interpolation method slightly to obtain 
(J7J) and then the slightly better lower bound OH]) on n. 
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